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ABSTRACT 

We study the halo mass function in the presence of a kurtosis type of primordial non- 
Gaussianity. The kurtosis corresponds to the trispectrum as defined in Fourier space. The 
primordial trispectrum is commonly characterized by two parameters, tnl and gNL- We focus 
on Tnl which is an important parameter to test the physics of multi-field inflation models. As 
applications of the derived non-Gaussian mass function, we consider the effects on the abun- 
dance of void structure, on early star formation, and on formation of the most massive objects 
at high redshift. We show that by comparing the effects of primordial non-Gaussianity on 
cluster abundance with that on void abundance, we can distinguish between the skewness and 
the kurtosis types of primordial non-Gaussianity. As for early star formation, we show that the 
kurtosis type of primordial non-Gaussianity seems on the average not to affect the reioniza- 
tion history of the Universe. However, at high redshifts (up to z ~ 20) such non-Gaussianity 
does somewhat affect the early stages of reionization. 
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1 INTRODUCTION 

The inflation paradigm has been well-known as a successful scenario for resolving several shortcomings of the standard Big Bang Model, 
in particular, the generation of primordial fluctuations which seed cosmic microwave background (CMB) fluctuations and structure forma- 
tion of the Universe. In the standard inflationary scenario, the primordial density fluctuations are generated from quantum fluctuations of 
a scalar field and they have almost Gaussian statistics. In recent years it has been realized that studying the non-Gaussianity of the pri- 
mordial density fluctuations can reveal valuable inf ormation about the dynamics of inflation (jlComatsu & Spergel 2001; Bartolo et al. lOoJ 
iBartolo. Matarrese & Riotto1l2010l : lKomatsij|20I0l) (and references therein). Thanks to significant progress in cosmological observations, 
most notably the CMB observations, we may expect that a meaningful measurement of this quantity will become observationally available 
in the near f uture and will thereby allow several inflation models to be tested. 

In Ref. jKomatsu & S pergef 200l|), the authors have introduced a simple new parameter which describes the deviation from Ga ussian- 
ity of the statistics of the primordial curvature fluctuations, the so-called non-linearity parameter /nl, defined as tSalopek & Bon^[l99d : 
iGangui et alJll994IVerde et al.l2000l) 

C(x) = Cg(x) + |/nl (CS(x) - {Cg(x)^)) + 0(CS(x)) , (1) 

where C, represents the primordial curvature fluctuations on a uniform energy density hypersurface and <^g denotes the Gaussian part. In the 
Probability Density Function (PDF) of the primordial fluctuations, the non-zero value of the non-linearity parameter /nl may generate a non- 
zero value of the skewness (3rd order moment), the kurtosis(4-th order moment) and so on. Obviously, the skewness can be parametrized by 
the leading term using /nl - However, the kurtosis can be affected not only by the /nl term but also by higher order terms, such as the Cg(''^) 
term in the above expression ([T). In general, one needs two parameters in order to characterize the kurtosis in the PDF. These parameters are 
normally called tnl and (?nl, where the first is usually (although not always) related to /nl and the second is the parameter that characterizes 
the third moment of C,. Such kind of non-linearity is the so-called local type of non-Gaussianity. Recently, other types of non-Gaussianity 
have been discussed in the literature, e.g., equilateral and orthogonal types. Theoretically, the local type of non-Gaussianity can be generated 
from the super-horizon non-linear dynamics of primordial curvature perturbations. On the other hand, the equilateral and orthogonal types 
of non-Gaussianity can be generated when one considers a scalar field which has a non-canonical kinetic term or the higher order derivative 
correction terms. In this paper, we focus on the local type of non-Gaussianity and consider the case where the equilateral and orthogonal 
types are negligible. 
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In the case where the primordial curvature fluctuations were generated from single field stochastic fluctuations (single-sourced case), 
i.e., the primordial curvature fluctuations can be expressed as Eq. l[l}, and tnl can be described only by /nl- But in general, e.g., if 
the pr imordial curvature fluctuation s were gene rated fr om multi-stochastic fluctuat ions then tnl and /nl have no universal relation any 
more jSuvama & Yamaguc hi 2008: Suvama et al.. 2010, , , Sugiyama. Komatsu & Futamase 201 1). Hence, it seems to be important to investi- 
gate the observational consequences of tnl independently of /nl- 

In this paper, we focus on the effects of the kurtosis type of primordial non-Gaussianity on the Large Scale Structure (LSS), 
in particular, on the halo mass function. There are many st udies of the effects of primordial non-Gaussianity on the LSS and also 
on the formulation of the non-Gaussian halo mass function /Matarrese. Verde & Jimenez 2000; Slosar et al.' 2008; 'Magg iore & Riottol 
I2OIOI ; IVerdell2010l ; iD'Amico et"al]|2010l iDe Simone. Maggiore & Riotto 2010; Wagner, Verde & Boubekeur.2010,) (and references therein ), 
which focus not only on /NL-type but also (/NL-type jDesiacgues & Seliakll20 1 d ; iMaggiore & Riottoll20 1 d ; I IChongchitnan & Siigi2010allbl ; 
lEngvist. Hotchkiss & Taanil3l20ld) . Here, we study the effects of kurtosis of the non-Gaussian primordial fluctuations whose non-linearity 
is parameterized by the two fr ee parameters, anl and tnl- Recently , a number of authors hav e studied non-Gaussian initial perturbations in 
two-field inflationary models jTseliakhovich. Hirata & Slosadl2010l ; ISmith & LoVerdell2010l) . In these papers, the authors have considered 
the effect of non-Gaussianity on the halo bias. Although this type of primordial non-Gaussianity is similar to the one considered here, we 
study the effects on the halo mass function. 

This paper is organized as follows. In the next section, we briefly review the kurtosis type of primordial non-Gaussianity considered 
here. In section |3] we formulate halo mass functions with primordial non-Gaussianity, based on the Press-Schechter theory and Edgeworth 
expansion. In section|4l we apply the non-Gaussian halo mass function to the formation of astrophysical objects. We consider three applica- 
tions; early star formation, the most massive object at high redshift and the abundance of voids. Section|5]provides a discussion and summary 
of our results. We adopt throughout the best fit cosmological parameters taken from WMAP 7-year data. 



2 TRISPECTRUM OF PRIMORDIAL NON-GAUSSIAN CURVATURE FLUCTUATIONS 

Here, we focus on the local-type non-Gaussianity. Following the notation commonly used, in the single-sourced case, up to the third order, 
the primordial curvature fluctuations can be expressed as 

C = Cg + |/nl (Cg - (Cg)) + ^5nlCS • (2) 
Based on this expression, the trispectrum of is given by 

(C(ki)C(k2)C(k3)C(k4)) = (27r)-'Tc(fci, /C2, fcg, fc*)^'" (ki + ka + ka + ki) , 

54 

rc(fci,fc2,fc3,fc4) = TNL (Pc('=i)'Pc(fc2)Pc(fci3) + llperms.) + —3NL(Pc(fci)-Pc(fc2)Pc(fc3)+ Sperms.) , (3) 

where kis = |ki + ksl and P(^{ki) is a power spectrum of C, given by (C(ki)C(k2)) = {2TvfP{ki)5^^'> (ki + kz). For the above definition 
of Tnl and the form of the non-linearity of the curvature perturbation (|2j, tnl can be written in terms of the non-linearity parameter /nl as 

-TNL = ^/nl ■ (4) 

This consistency relation is satisfied only in the case where the primordial curvature fluctuations can be described by Eq. l[2ll, namely, the 
primordial curvature fluctuations are sourced only fr om the quantum fluctuations of a single scalar field, e.g., curvaton jEnavist & SlothI 
l2002l ; lLvth & Wandsll2002l ; lMoroi & Takahashill200lh . 

However, if there are mu l tiple sources of the primordial curvature fluctuations, th en the above consistency relation is not satis- 
fied l lLanglois & Vemizzill2004l; llchikawa et"al]|2008l ; lHuan3l2009l ; Isymes & Choi Il20ld)^ In general, it has been known that there exists 
an inequality between the local type n on-linearity parameters tnl and /nl given bv dSuvama & Yamaguchi| [200a; ISuvama et alj [201ol; 
ISugivama. Komatsu & Futamasdl2oTlh 

TNL > - -/nl • (5) 



For example, let us consider the local-type non-Gaussianity given by 

C = (l>G+ |/nl (<^g - {(Pg)) + tK-L(t)Gi>G , (6) 

where (I>g and are Gaussian fluctuations with {(f>GtpG) ~ and tNL is a non-linearity parameter, which represents the non-linear coupling 
between (pa and xpc in (g ■ At leading order, the power spectrum of is given by that of the Gaussian part (pa as 

(C(k)C(k')) = (0G(k)0G(k')) = (2^)=^P4fc)5'''(k + k') , (7) 
and the bispectrum is given only by /nl as 

{C(ki)C(k2)C(k3)) = (27r)='|/NL(P^(fci)P<^(fc2) + 2perms.)<5(')(ki + ka + kg) , (8) 
because of {(pG^c) ~ 0. 



Modification of the halo mass function by kurtosis associated with primordial non-Gaussianity 3 

In the single-source case which corresponds to the case of tNL = 0, as mentioned above, the trispectrum can be also parameterized only 
by /nl- However, for the above type of curvature fluctuations the trispectrum is given by 

Tc(fci,fc2,fc3,fc4) = f|/NL^ (P^(fci)P^(fc2)P<^(fci3) + llperms.) 



5-' 

P</,(fcl)P'0(fc2)Pv'(fcl3) +f'0(fcl)^'0(fc2)P^(fcl4) +P^(fcl)P0(fc3)P^(fcl2) 
+P4.{k2)P4.{k3)Pi,{kl2) + P^{k2)P4.{k3)Pi,{k2i) + P0(fc2)P</,(fc4)P^(fcl2) 

+P^{k2)P^{k4)P^,{k2-i) + P4.{k3)P4,{ki)P^,{kri) + P0(A:3)P<^(fe4)P^(fc23) ) , (9) 



where fcis = |ki + k3|. We assume that the power spectra of random Gaussian fields (f)G and ipa have only weak scale-dependence, that is, 
the power spectra are respectively given by 

'^^('^^-l^^HW ' '''' 

where ka is a pivot scale and — 1| <C 1 and — 1| <C 1. In such a case, we can rewrite the power spectrum of ?/)g as 

P^[k) ~ oiP^{k) , a = A^/A^ , (II) 

and then using the ratio of the amplitudes a, the expression for the trispectrum can be reduced to 



{C(ki)C(k2)C(k3)C(k4)) ^ (27r)^5(^)(ki k2 + kg + k4) ( ^/^L + q4l ) (P,/. (A:i)P4fc2 )P0 (few) + llperms.) 
From the above equation and Eq. (|3}, we easily find that the non-linearity parameter tnl is 



(12) 



'^^ -c2 , ,2 ^ 36 2 

TNL = ^/nL + «InL ^ ^^NL • (13) 

Hence, in the following discussion, we consider tnl independently of /nl- 



3 NON-GAUSSIAN MASS FUNCTION INDUCED FROM PRIMORDIAL NON-GAUSSIANITY 

In the previous section, we have shown that there is a strong theoretical motivation for considering tnl to be independent of /nl. The 
parameter tnl characterizes the amplitude of the trispectrum of primordial curvature fluctuations as well as ^nl- Here, we briefly review the 
formula for the halo mass function with not only the non-zero primordial bispectrum but also the non-zero primordial trispectrum, based on 
Press-Schechter theory. 



3.1 Probability Density Function of the smoothed density field with primordial non-Gaussianity 

The matter density linear fluctuations in Fourier space at redshift z, S(k, z), are given by the primordial curvature perturbation on a uniform 
energy density hypersurface ("(k) as 

SCk,z)=M{k)D{z)Cik), (14) 

Mik)^l-^^T{k), (15) 

O iZmO J^o 

where flmo is the present density parameter for total non-relativistic matter. Ho is the Hubble constant, D{z) is a linear growth function and 
T{k) is a transfer function. Using these expressions, we can obtain the linear matter power spectrum as 

{5{k,z)S{k\z)) = (27r)=*Pi(fc,z),5<^'(k + k') , (16) 

Ps{k, z) = ^M{kfD{zfV4,{k) , (17) 

where V^{k) = k^ P^{k) /{2it^). Following the standard procedure, let us define the smoothed density fluctuation on a given length scale, 
R, as 

(2^ 



^« = / 7:rT3Wfl(fc)^(k,z) , (18) 
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where Wi?,(fc) is the Fourier transform of a spherical top-hat window function given by 



(19) 



In order to take into account primordial non-Gaussianity in the smoothed density fluctuations, let us consider the PDF of Sr, F{Sfi)d5u. 
The n-th central moment for F{Sfi)d5u is defined as 



(Sn) 



5lF{5R)d5R 



and each reduced p-th cumulant can be defined as 



S^{R) 



(20) 



(21) 



where a subscript c denotes the connected part of p-point function given by 

{5r)c = , {5'r)c = {5\) = cj\ , 

(4)c = (4) , (4)c = (4) - 3{5j,)l 



etc.. 



(22) 

with zero mean density field. Here, cr^, S3{R) and S4{R) are the variance, the skewness and the kurtosis, respectively. Let us consider a 
non-Gaussian PDF of matter density fluctuations, based on the concept of the Edgeworth expansion. Here, we co nsider the expansion of 
the PDF of the dens ity field F{i')du with v = Sr/ur, in terms of the derivatives of the Gaussian PDF, Fg{v), as jjuszkiewcz et al.ll 19951 : 
ILo Verde et all2008i) 



F{v)dv — dv 
with 



Co 



(27r)-'''2exp(-i/V2) 



F^^\u) 



d 



-, — Faiu) = 



irHmiu)FGiu) 
where Hm^v) is the Hermite polynomials; 



3u . 



(23) 

(24) 
(25) 

(26) 



From the above relation between the derivatives of the Gaussian PDF and Hermite polynomials, we can regard the expression ([23} as a 
non-Gaussian PDF expanded in terms of the Hermite polynomials. Since the Hermite polynomials satisfy orthogonal relations; 



H,4u)H„{u)FG{u)du = 
we can evaluate the coefficients as 







, if m 7^ n . 
, if m = n . 



Hm{v)F{v)dv . 



(27) 



(28) 



Then, we can obtain the expressions for the coefficients, Cm, in terms of the reduced cumulants (variance, skewness, kurtosis and so on) as 

Co = 1 , Ci = C2 = , C3 = —S3{R)(JR , C4 = Sa{R)(j1i , 

C5 = -S5{R)a% , C6 = lOSaiRfal + S6{R)a% , • • • , (29) 
and, as a result, the non-Gaussian PDF of the density field, F{u)dv, can be obtained as 



F(v)dv = 



dv 



exp ( 



V2) 



S;{R)(jR 
6 



H,{v) + - 



S3{R)(jR 
6 



24 



24 



24 



(30) 



up to the third order terms in Sz{R) and S4,{R) and neglect the contributions of the higher order cumulants; Sn{R) {n ^ 5). This derivation 
of the non-Gaussian PDF is based on the so-called Edgeworth expansion. Of course, the non-zero non-linearity parameters /nl, tnl and 
(?NL also generate non-zero higher order cumulants; Sn{R) {n ^ 5). However, as far as considering the non-Gaussian curvature fluctuations 
given by Eq. l|6) and current observational constraint s on the non-linearity parameters ( Komatsu et al.ll201 H ; iFergusson. Regan & Shellard 
terms Sn{R) {n ^ 5) are greatly suppressed dEnqvist. Hotchkiss & Taanilall2010h . Hence, the assumption of neglecting the higher 
order cumulants seems to be reasonable. 
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3.2 Halo mass function with non-Gaussian corrections 



Let us consider the halo mass function with non-Gaussian PDF of the smoothed density field as given in the previous subsection. Based 
on the spirit of the Press-Schechter formula, the halo mass function which gives the number density of coUapsed structures (halos) with 
the mass between A/(= 47rp_R^/3 with p is the background matter density) and M + dM at a redshift z, (dn(M, z) / dM)dM is given 
by l^'Amico et al. 201(D 



dn 
dM 



(M, z)dM = 



'dM 



= ~dM 



2p d 
MdM 

2 p 



duF(u) 



-yyexp 

TV Al 



' 2 



d In an 



+ 



6 

Si{R)a%, 



dM 
1 f SsjR) 



+ 



2 V 



6 



+ 



24 

d ( S3{R)CJR 



dM\ 6 



Ha{v, 



1 ( Si{R)al 

2 V 24 



6 V 6 

1 (SA{R)al 



6 



1 d 
2dM 



dM 



24 



S3iR)aB. 
6 

Si{R)a% 



24 

1 d ( Sz{R)aR 



24 



6dM V 



6 

g4(J?)a| 
24 



ifii(i^c)|+-- - (31) 

where = Sc/o-r and (5c denotes the threshold for collapse which is originally given by 5c ~ 1.69. However, in Ref. l lGrossi et al]|2009l) . 
the authors have suggested that using the correction 5c — >■ (5c with q — 0.75 puts the analytic predictions in good agreement with the 
numerical simulations. This is due to the more realistic case of ellipsoidal collapse. Hence 5c — 1.69 x is often referred to as the critical 
density of ellipsoidal collapse. Here we adopt this corrected density threshold Sc — 1.69 x \/0.75. In the following calculations, we use the 
above formula of the non-Gaussian mass functions up to the third order in terms of 5*3 and S4. 
For a Gaussian probability distribution, the mass function is given by 



dno 
'dM 



(M, z)dM 



[2 


P 




V TT 


M 


2 



dM 



VcdM 



(32) 



(34) 



and we define the ratio between the non-Gaussian mass function and the Gaussian one a|3 

^ , dn{M,z)ldM 

Let us focus on the redshift dependence of the above expression. From the definition of the reduced cumulants M\\ and the fact that the 
redshift dependence of the density field is given by 5r oc D{z), we can easily find that a^^'^ Sp{R) has no redshift-dependence. Hence, any 
remaining redshift dependence comes only from the term Here, following the literature, the redshift-dependence can be carried by 5c as 
5c — > 5c{z) oc D{z)^^ and then the variance aR has no redshift-dependence. In the following discussion, we change the subscript Rto M 
because R and M have a one-to-one correspondence through the equation M — AtvR'^P/S. 



3.2.1 Variance, skewness and kurtosis 

Let us consider the concrete expressions of the variance, skewness and kurtosis of the primordial curvature perturbations whose power-, bl- 
and tri-spectra are given by Eqs. Q, (|8) and (|3), respectively. The variance is given by 

f dk 

^1= / ^WUk)Mikfr^{k) , (35) 



1 In Refs. iinqvist. Hotchkiss & TaanilalioToljMatarrese. Verde & Jimenezl200(]| : IVerde et al.l200ll) . the authors introduced the MVJ convention for defining 
the ratio given by 



cxp 



3S3{R)CTR , 4S4(i?)a2 



24 



X 



S:i + 



yc ( Si{R)aR\ , (dlncTRY^ d /SsiKjoR 



S3 



6 J ^ dM J dM 6 



54- 



[ S4{RWr\ /dlnaR\-^ d ( Sa{FI)^1 



54 



12 



dM 



dM \ 24 



<53 = 



1/2 



(33) 



which is not based on the Edgeworth expansion. In our calculation, we have also checked the consistency between the above MVI expression and Eq. i34\ . 
This issue is discussed in Appendix IaI 
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z II ki 



X 




Figure 1. the thi'ee vectors, ki, k2 and kg in the trispectram. 



the skewness is dChongchitnan & SilldbOlOiJ) 

S3{R) = 



dfj,i 



■WR{ki2)M{ki2) 



P4.(ki 



P^{k2 



(36) 



where ki2 
by 

SUR) = 
SUR) = 



^/fcj + fcj + 2fcifc2/ii2 and jj.12 — cos 612, and the kurtosis which is proportional to the non-linearity parameter tnl is given 



^W/fl(fci)A1(fci)P4fcl) f ^WR{k2)M{k2)V^{k2) 
fel J fe2 

-P,^>(A;i23 



Wfl(A;3)X(A:3)P<^(fc3) 



X<^ P^(fcl2) 



-P0(fc23) 



1 



+ 



+ 



P4,{k2 

1 



1 + 



P0(fc2) P0(fc3) 



P</,(fcl23) 



P,^,(A:i 



+ P4.{k3 



1 1 

P^(A:3) ^ P^(A:i 



1 + 



P</,(fci 



P,#,(fc2 



(37) 



Here, we have fixed the three vectors, ki, k2 and ks that appear in the expression of the trispectrum, as shown in Fig.[T] Hence, using the 
angular variables, 612, 613 and ipis, we have 



fel2 = \/kf + kl+ 2fclfc2/il2 



k23 = y'fc| + fc| + 2fc2fc3 (1 - Ml2)(l - fJ-ls) COS Ifiis + /il2Ml3^ , 



fcl3 = 



fel + fe| + 2felfe3/X13 



and 



fel23 



fcf + fci + fe| + 2felfc2A112 + 2fel 



fesMlS + 2fe2fc3 (1 — Ml2)(l ~ M13) COS V?13 + ^12/il3^ , 



(38) 



(39) 
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where ^tj = cos 6ij . In order to calculate the skewness more easily, let us consider the squeezed limit in momentum space, e.g., fci <c; ^2 — 
fcs . In this limit, the equation for the skewness ( I36t can be reduced to 

S3 :^2a% f ^WR{ki)M{ki)p4,{ki) , (40) 

and by considering other limiting cases, i.e., ^2 ^> and kz 0, we obtain 

S3 ~ 6ai j ^WR{k)M{k)V4,{k) . (41) 
Based on the above approximate expression, we find a simple formula; 

(^RSziR) = 4.3 X 10~Vnl X cr^-'^ (IO'^/i^'Mq < M < 2 x 10'^/i"^Mq) . (42) 

This result seems to be close to those given in Refs. jPe Simone. Maggiore & Riottol20Iol ; lEnavist. Hotchkiss & TaaniIal201oH^ . Hence, we 
adopt the above expression in the following discussion. In a similar way, from the expression of the kurtosis ( 137) , we can easily find that the 
kurtosis induced from the non-linearity parameter tnl becomes largest in the limit of ki Q(i — 1, 2, 3, 4) or kij — >■ 0{i 7^ j = 1, 2, 3, 4) 
(local type). Then, we have an approximate expression 

/Hk 
— WR{k)M{k)V4k) X SsiR) + l2A^a% . 

(43) 

On the other ha nd, in the squeezed limit kj — >■ 0(^ — 1, 2, 3, 4), the kurtosis whic h is proportional to the non-linearity parameter (?nl can be 
also reduced to jChongchitnan & Silkl2010al : lEnqvist. Hotchkiss & Taariii3l201(]h 



•^4 — TTF^Ze"''' 



54 5nl ?,g 
257^^^ ' 



'R 



— WR{k)M{k)V^{k) X S-s{R) . (44) 

From these approximate expressions, we respectively obtain simple formulae for the kurtosis in the form 

a%Sl{R) = 1.9 X IO^Vnl X cr^^^ (IO'^/i^'Mq < M < 2 x IO'^^/i^'Mq) , 

cjrSUR) = 9.4 X 10"**gNL x a^^^ (10'^/i~'Mq < Af < 2 x 10'^fe"'Mo) . (45) 

The result for 5| (i?) also is close to that obtained in Ref. jEnqvist. Hotchkiss & Taanilall2010lFl Hence in the following discussion , we also 
adopt the above expressions for the kurtosis as well as that for the skewness. 



3.2.2 Dijference between the Gaussian and the non-Gaussian mass functions 

Based on the above calculations for the variance, o"^, the skewness, S3, and also the kurtosis, S4, the mass function can now be calculated. 
In the following discussion, we take values of the non-linearity parame ters as /nl ~ 10 0, tnl = 10® and ^nl ~ 0. This value of tnl 
may be inconsistent with the observational constraint obtained by Ref. (Smidt et al. '2010') as -0.6 < tnl/IO" < 3.3 at 95% confidence 
level. However, there might be a caveat since in Ref. ( iFergussonTR egan & Shellard 2010), the authors have claimed that the approach in 
Smidt et al. does not directly subtract the effect of anisotropic noise and other systematic effects which are important in obtaining an accurate 
and optimized result. Nonetheless, in order to emphasize the differences between the Gaussian mass functions and the non-Gaussian mass 
functions with the non-zero /nl and the non-zero tnl cases, we take the above values. 

In Fig.|2l we show that the mass function in the mass range between 5.0 x IQ^'^H'^Mq and 2.0 x lO^'"^ H'^Mq at the redshift z = 0. The 
red thin line shows the mass function with the Gaussian density fluctuations given by Eq. ( |32| (. The blue dashed and green thick lines show 
the non-Gaussian mass function given by Eq. l l31b in the cases with /nl = 100 and tnl = Snl = and tnl = 10® and /nl ~ 3nl ~ 0, 
respectively. From this figure, it is rather difficult to see the differences between the Gaussian and the non-Gaussian mass functions. In Fig. [3] 
we show the ratios between the Gaussian and the non-Gaussian mass functions defined by Eq. ( |34t . The red dashed line shows TZng{M, 0) 
with /nl ~ 100 and tnl ~ .gNL ~ and the blue solid line shows that with /nl = .gNL ~ and tnl = 10®. The magenta dotted line i s 
for the case with tnl ~ ffNL ~ and /nl ~ 30 which is corresponding to the mean value of the current WMAP data jKomatsu et alj2oTll ). 
The black dashe d-dotted line is for the case with /nl ~ (Jnl ~ and tnl ~ W'^ which is consistent with the maximum allowed value 
obtained by Ref. dSmidt et alj2010l) . From this figure, we infer that for both types of primordial non-Gaussianity, i.e., positive skewness and 
kurtosis, the mass functions can be systematically enhanced for more massive objects, as compared with the Gaussian case. The enhancement 
of the mass functions depends on the values of tnl and .gNL. We find that for the cases with /nl ~ 30 and tnl = 10** 7?.ng are respectively 
1.06 and 1.01 for M = 2 x 10^^ H'^Mq. Hence, in both cases the effects of the primordial non-Gaussianity on the mass functions seem to 



As mentioned in Ref. iEngvist. Hotchkiss & Taanilj|2010l) . this result is different from that in the published version of Ref iChongchitnan & Silkll2010al) . 
However, the authors in Ref. iChongchitnai7^Sil32010al) have corrected the result in the arXiv version and their new derivation is now close to our result 1 1421 . 

In addition to the expression derived for the skewness, the result for kurtosis in Ref. tChongchitnan & Silk.2010ab has been corrected in the arXiv version 
and this is also close to our result i45\ . 
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Figure 2. The mass function with the mass range between 4.5 X lO^'^h"^ Mq and 2.0 X lQ^^h~^ Mq at the redshift z = 0. The red thin hne shows the 
mass function with the Gaussian density fluctuations given by Eq. )32t . The blue dashed and green thick lines show the non-Gaussian mass function given by 
Eq. i3i\ in the case with /nl = 100 and tkl = 3nl = and the case with tnl = 10^ and /nl = 9NL = 0, respectively. 




12.5 13.0 13.5 14.0 14.5 15.0 

log.o {M/[h-'Me]) 



Figure 3. The ratio between the Gaussian and the non-Gaussian mass functions. The red dashed line shows T^ng {M, 0) with /nl = 100 and ttml = 9NL = 
and the blue solid line shows that with /nl = 9nl = and tnl = 10^. The magent a dotted line is for the case with rj^L = 9nl = and /nl = 30 
which is corresponding to the mean value of the cuiTent WMAP data iKomatsu et alj20Tl[). The black dash ed-dotted line is for the case with /nl = 9NL = 
and T>jL = 10* which is consistent with the maximum allowed value obtained by Ref. Ismidt et alj2010t) . 

be too small to detect. We also find that the enhancement of the non-Gaussian mass function with the non-zero kurtosis type of primordial 
non-Gaussianity, i.e., non-zero tnl, depends more strongly on the mass of the collapsed objects than the case with the non-zero skewness 
type of primordial non-Gaussianity. This is because in the expression for the non-Gaussian mass function ( |3U . the (5c /o-_h -dependence of the 
term related with the kurtosis S4 is stronger than that of the term related with the skewness Sa, namely, S'4-term oc {5c/<jr)^ and Ss-term 
oc {Sc/ctrY- As the collapsed objects become more massive, the variance or becomes smaller and hence 5c/<yR becomes larger. Thus, if we 
would detect the enhancement of the mass function for massive collapsed objects and find its scale-dependence, then we might distinguish 
the kurtosis type of primordial non-Gaussianity from the skewness type. In Fig. |4l we show the redshift dependence of the ratio between 
the non-Gaussian mass function and the Gaussian mass function as we change the value of tnl. Here we have fixed the mass of the halo 
as M = 10"/i"^A/o. The solid line is for the case with tnl = 10^, the dashed hne for 2.0 x 10^, dotted line for 5.0 x 10^ and the 
dashed-dotted line for 10®. From this figure, we find that at higher redshift the enhancement of the mass function for massive collapsed 
objects increases. This is because the critical density 5c{z) = 5c/D{z) becomes much larger at larger redshifts due to the smaller linear 
growth function D{z). Hence, in order to observationally test the kurtosis type of primordial non-Gaussianity it will be useful to observe 
high-redshift rare objects. 
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Figure 4. The ratio between the non-Gaussian mass function and the Gaussian mass function v.s. the redshift (0 < z < 5) with changing the value of tnl- 
Here we have fixed the mass of halo as M = 10^'*/i~^Mq. The solid line is for the case with tnl = 10^, dashed line for 2.0 X 10^, dotted line for 5.0 X 10^ 
and dot-dashed line for 10®. 



4 APPLICATIONS 

In this section, we consider applications of the mass function with both skewness and kurtosis types of primordial non-Gaussianity. Here, we 
also take values of the non-linearity parameters to be /nl ~ 100, tnl = 10^ and ^nl ~ 0. 



4.1 Early Star Formation 

Let us first investigate the effect of primordial non-Gaussianity on the epoch of reionization. As is well-known, in ord er to understand the 
mechanism of reionization, it is important to estimate the number of photon s from Population III stars. Following Refs. dSomerville & Liviol 
I2OO3I : ISomerville. Bullock & Liviol l2003l : |s ugivama, Zaroubi & Si"iMl2004h . the global star-formation-rate density denoted by p, can be 
written as 

p. = e*p6^i?h(Mvir > M > Merit, t) . (46) 

Here, pj, is the background baryon number density and e, denotes the star-formation efficiency usually taken to be 0.002 for 2OOM0 Pop 
III stars and 0.001 for WOMq. Fh{Mvii > M > Merit, t) represen ts the fraction of the total mass in collapsed ob jects (halos) with masses 
greater than the minimum collapse mass scale Merit = lO^/i^^M© jVoshida et al.l2003l ; lFuller & Couchmanll2000l) and lower than the virial 

mass Mvir = M{T^ir = WK). 

The relation between the mass and the virial temperature is given by jSarkana & Loeb ll200ll : lYoshida et ai]|2003l) 

where p is the mean molecula r weight, and Ac{z) is the final overdensity relative to the critical density, which is given by a fitting for- 
mula tervan & Norman Ill998h 

Ac = IStt' + 82 {n^{z) - 1) - 39 (n„(2) - 1)' , (48) 
where Qrn{z) is the density parameter of matter at redshift z; 

Assuming that the photon number production-rate per Mq from Pop. Ill stars is N^, = 1.6 x 10**s~^Mg ^ and that the life time of Pop. Ill 
star is rm = 3.0 x lO^yr, we can obtain the total production rate of ionizing photons at time t as 

^(t) = e.pftiV^ {Fn{t) - F,{t - r„i)) , (50) 
hence the cumulative number of photons per H atom is 

^{z) ~ pmpe,iV^Fh(Mvir > M > A/crit, z)rui , (51) 
with the proton mass nip and the hydrogen number density uh- In the above expression, Fh{Mvir > M > Merit, z) is given by Press- 
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Figure 5. Cumulative photon number per H atom as a function of tfie redsliift for 8 < z < 18. The black solid line shows n^/riH (z) for the case with the 
Gaussian fluctuations. The blue dashed line is for the case with the non-Gaussian fluctuations; /nl = 100 and tnl = Snl = 0. The red dot-dashed line is 
for the case with tjml = 1-0 x 10*^ and /nl = 9nl = 0. The thin black dashed line corresponds to n-y/riH as a guide of the complete reionization. 
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Figure 6. We plot the ratio between n^(z) /nu in the pure Gaussian primordial fluctuation case and that in the non-Gaussian case for 5 < z < 20. The blue 
dashed line is for the case with /nl = 100 and tnl = 9nl = and the red solid line for the case with t^l = 1-0 x 10^ and /nl = Snl = 0. 



Schechter theory as 

Fh{M^i. >M> Merit, 2) = - M^(M, z)dM . (52) 

Substituting our expression (|3D for the non-Gaussian mass function into the above equation, we can estimate the effect of primordial 
non-Gaussianity on the number of photons emitted from Population III stars, which is one of the most important quantities during the epoch 
of reionization. In Fig.|5] we show the cumulative photon number per H atom given by Eq. l |5Ib as a function of the redshift for 8 < z < 18. 
The black solid line shows n-f/nH{z) for the case with the Gaussian fluctuations, the blue dashed line is for the case with the non-Gaussian 
fluctuations; /nl ~ 100 and tnl = .gNL = and the red dot-dashed line for the case with tnl = 1.0 x 10® and /nl = <7nl = 0. The thin 
black dashed line coiTesponds to n-y/nn = 10 as a guide of the complete reionization on average. From this figure, we find that primordial 
non-Gaussianity seems not to affect the reionization history of the Universe on average which is characterized by the value of n-y/nn ~ 10. 
However, at higher redshift the effect of primordial non-Gaussianity seems to be significant on the cumulative photon number density. We 
evaluate this effect in Fig. (6) where we plot the ratio between n-y/nniz) in the pure Gaussian primordial fluctuation case and that in the 
non-Gaussian case. The blue dashed line is for the case with /nl = 100 and tnl = Qnl ~ and the red solid line for the case with 
T"NL = 1.0 X 10® and /nl ~ 3nl = 0. From this figure, we find that compared to the Gaussian case the cumulative number of photons in 
the non-Gaussian case is larger at higher redshifts both in the non-zero S3 and the non-zero S4 cases. Moreover, as we have mentioned in the 
previous section the kurtosis type of primordial non-Gaussianity affects the enhancement of the photon number density more significantly 
at high redshift. That is, there seems to be the possibility of dramatically changing the history of the early stage of reionization due to the 
kurtosis type of primordial non-Gaussianity even for values in the range of the current limits obtained from CMB observations. Of course, 
the above rough estimate is not precise enough to enable us to estimate the exact cumulative number of the ionizing photons. However, we 
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consider here that, in view of the completely ad hoc nature of the amount of non-Gaussianity due to the absence of a compelling inflationary 
model, it suffices for us to focus on the deviation of the photon number based on the non-Gaussian mass function from that based on the 
Gaussian mass function. 



4.2 High-Redshift Massive Clusters 

Recently, the authors in Ref. jjee et al.l2009l;lRosati et al.l2009l) have presented a weak lensing analysis of the galaxy cluster XMMU J2235.3- 
2557 which has a high redshift z ~ 1.4 and whose mass is M324 = (6.4 ± 1.2) X lO^MeB In ACDM model the formation of such a 
massive clus ter at this redshift would be a rare event (at least 3cr). 

In Ref. jCavon. Gordon & Sil]dl201ol) . the authors have considered the effects of primordial non-Gaussianity parametrized by the non- 
linearity parameter /nl which they found to be /nl ~ 449 ± 286 at wave number of about 0.4Mpc^^ in order to explain the existence 
of such a massive cluster at high redshift. Considering scale-invariant /nl, this result contradicts the current CMB observational constraint 
JwL< 100. Therefore, the authors remarked that one would need to invoke scale-dependent /nl. In Ref. tengvist. Hotchkiss & Taanil3 
I2OIOI) , the authors have considered non-zero jjnl case and found that (;nl = 0(10*^) could explain the existence of high redshift massive 
clusters. 

Here, instead of considering the scale-dependence of /nl or qnl, let us consider the effect of the kurtosis induced from the TNL-type 
primordial non-Gaussianity on the formation of massive clu sters. Of course, for a more detailed analysis we need to calculate the probability 
of the massive clusters under the procedure done in Ref. dCavon, Gordon &ll ilk 2010). However, in order to give a naive estimation of 
the value of tnl which can explain the existence of the massive cluster XMMU J2235. 3-2557, we investigate the value of tnl which 
gives the same value as does the non-Gaussian mass function, namely, TZng defined as Eq. ( I34t . including the effect of kurtosis S4 on the 
corresponding scale at the corresponding redshift by including the effect of /nl, i.e., skewness. Here, we adopt Mxmmu ~ 6.4 x lO^^M© 
and 2XMMU ~ 1-4 as the mass and the redshift of the massive cluster XMMU J2235. 3-255, respectively. For the value of /nl, the best fit 
value derived in Ref. dCavon. Gordon & Silkll2010l) is adopted. For these parameters, we also find that this value can be realized in the case 
with /nl ~ .9NL = and tnl = 1.7 x 10^. As we have mentioned in Sec. 13.2.21 this value may be ruled out by the result obtained by 
Ref. I Smidt et al]|2010h . Hence, if we believe this constraint, we need to consider the possibility such as scale-dependent tnl. 



4.3 Abundance of voids 

As an other example, we study the void abundance with primordial non-Gaussian corrections. In Ref. jKamionkowski. Verde & Jimenez! 
l2009l) . the authors showed that the void distribution function can be derived in the same way as the halo mass function using Press-Schechter 
theory. This is done by replacing the critical "overdensity" parameter, Sc, with the negative "underdensity" parameter, The precise value 
of Sv depends on the definition of a void. For example, if the voids are regio ns having a density half of p , then we can estimate the critical 
value of underdensity as Sy —0.7 jKamionkowski. Verd e & Jimenez 2009). There are also several numerical studies about the value of 5^ 
which suggest 5„ ^ -0.8 jShandarin et al.ll2005HPark & Lee 2007; Colberg et al. 20oH) . 

In any case, based on Press-Schech ter theory, the abundance of voids which have radius between R and R + dR is given 
by jKamionkowski. Verde & Jimenezll2009t) 
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Up to the third order in terms of S3 and ^4, the void abundance with primordial non-Gaussian corrections is also given by 
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^ The halo is defined as a spherical overdense region whose density is 324 times the mean matter density of the Universe. 




Figure 7. The ratio between the void abundance distribution function with non-zero primordial non-Gaussian correction and the Gaussian distribution function 
at redshifts z = 0.0, 0.75, 1.5 given by Eq. )56t . The red lines are for the case with /nl = 100 and tnl = (/nl = 0, the blue thick lines for the case with 
TjjL = 10*' and /nl = 9NL = and the black thin lines for the case with tnl = 10'* and /nl = 9NL = 0, The solid lines are for the case with z = 0.0, 
the dot-dashed lines for the case with z = 0.75 and the dashed lines for the case with z = 1.5. Here we took = —0.7. 



Following the previous section, we define the ratio between the void abundance with the pure Gaussian PDF and that with the primordial 
non-Gaussian corrections as 

oM ^ dn^°"'{R)/dR 
~ dnl?'''{R)/dR ■ 

(56) 

In Fig. [T] we show 7?.ng^ for the cases with /nl ~ 100, tkl = (/nl ~ (red lines), tnl ~ 10^, /nl = Qnl = (blue thick lines) and 
tnl = 10*, /nl ~ Pnl ~ (black thin lines). We adopt 5y = —0.7. We also show the ratio with changing the redshift; the solid lines for 
the case with z — 0.0, the dot-dashed lines for the case with z = 0.75 and the dashed lines for the case with z = 1.5. From this figure, we 
conclude that the non-Gaussian void abundance with non-zero tnl becomes larger than the Gaussian one on relatively larger scales whereas 
that with the non-zero /nl becomes smaller. On the other hand, as seen in the previous section, the halo abundance becomes larger not 
only with non-zero tnl but also with non-zero /nl in relatively more massive objects. Hence, from this discussion, we confirm that the 
non-Gaussian effects on both t he halo and the void abundan ces allow as to distinguish the large kurtosis, i.e., large tnl, case from the large 
skewness, i.e., large /nl case dChongchitnan & SilklbOlOj) . 



5 SUMMARY AND DISCUSSION 

It has recently become clear that cosmological large-scale structure and CMB observations could provide stringent constraints on the PDF 
of primordial adiabatic curvature fluctuations. In particular, the high order moments of the PDF, such as its skewness and kurtosis, can give 
unique insights into the dynamics and conditions of the inflationary phase in the early Universe. 

In this paper, we have investigated the effects of the TNL-type of primordial non-Gaussianity on the halo mass function. In particular, 
we have obtained a formula for the halo mass function with the non-Gaussian coirections coming from the kurtosis induced by the non- 
zero Tnl- We find that the deviations of the non-Gaussian mass function from the Gaussian one become larger for larger mass objects 
(M > lO"/i"^M0 for 2 ~ 0) as well as at higher redshifts {z > 1 for M ^ lO^/i'^M©) in the case with tnl = O(10'^). Such 
features are quite similar to those obtained from skewness-driven non-Gaussian corrections that are induced by the /NL-type of primordial 
non-Gaussianity. 

As examples of applications of our formulae, we have considered the effects on early star formation, formation of the most massive 
objects at high redshift, and the abundance of voids. 

For early star formation, we applied our formula for the non- Gaussian halo mass function in order to estimate of the redshift-dependence 
of the cumulative number of photons emitted from population III stars, a crucial quantity in considerations of the reionization history of the 
Universe. We found that primordial non-Gaussianity does not affect the reionization history of the Universe on the average, but at high 
redshift (z ~ 20), namely the earliest stages of reionization, it is effective. 



Modification of the halo mass function by kurtosis associated with primordial non-Gaussianity 13 



We have also obtained an estimate of the value of tnl needed to naturally explain the existence of the galaxy cluster XMMU J2235.3- 
2557, namely tnl = 1.7 x 10^. Hence, in li ght of the result of Smidt et al., we might need to consider a possibility such as scale-dependent 
TNL in the case with non-zero /nl. In Ref. faovle. Jimenez & Verddl20ld) . the authors have investigated 15 high-mass and high-redshift 
galaxy clusters and found that such objects are extremely rare in the standard ACDM model with Gaussian primordial fluctuations. They 
derived a constraint on /nl in order to explain the mere existence of th ese objects as /nl > 475 at 95% c onfidence level, with the other 
cosmological pa rameters fixed to best fit values of WMAP data. In Ref, tenqvist. Hotchkiss & Taanil3l2010h . the authors have extended the 
analysis of Ref. l lHovle. Jimenez & VerdelbOld) to the case with non-zero (?nl. It should clearly be of interest to derive a constraint on tnl 
for these observed high-mas s and high-redshift galaxy clusters. We will address this in future wor k. 

As mentioned in Refs. ( iKamionkowski. Verde & Jimenezl200Sl ; IChongchitnan & Silkl2010ah . the non-Gaussian correction coming from 
skewness reduces the abundance of voids on large scales when the non-linearity parameter /nl is positive in contrast to the fact that positive 
/nl enhances the number of more massive halo objects. On the other hand, the non-Gaussian correction coming from kurtosis enhances 
not only the numbers of more massive halo objects but also the abundances of voids on large scales. Hence, if one could also measure the 
void abundance as well as the halo mass function more precisely, one could potentially distinguish between the /nl and the TNL-types of 

primordial non-Gaussianity. 

NOTE; During the time that we were preparing this manuscript, Ref. jLo Verde & SmittJ 1201 ll) appeared on the arXiv. In 
Ref. jLo Verde & Smithll201ll) . they considered the same type of primordial non-Gaussianity as in our study and obtained a useful ana- 
lytic formula for the halo mass function with the kurtosis type primordial non-Gaussianity using N-body simulations. We find that our 
formula ( 131) is in reasonably good agreement with their formula as far as the behavior of the halo mass function with the kurtosis type of 
primordial non-Gaussianity. 
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APPENDIX A: MVJ EXPRESSION 

In Ref. jMatarrese. Verde & Jimene3l2000h . the authors have given a formula for the ratio between the non-Gaussian mass function and the 



Gaussian mass function as 
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which was not derived based on the Edgeworth expansion as mentioned in section [5] In section |4l we have discussed some applications of 
the non-Gaussian halo mass function. 

As for the discussion in subsection l4. 1 l about early star formation, the redshift-dependence of the critical value of the cumulative photon 
number per H atom (n^/nH = 10) is not so sensitive to primordial non-Gaussianity. In order to estimate more precisely how n~,/nH at 
high redshift is enhanced due to primordial non-Gaussianity, we should check which formula better describes the effect of primordial non- 
Gaussianity on the halo mass function. This is a future issue. In subsection 14.31 we have discussed the void abundance and noted that the 
kurtosis type of primordial non-Gaussianity can enhance the abundance of the large voids as opposed to the skewness type of primordial 
non-Gaussianity. This is just qualitative discussion. 

On the other hand, the discussion in subsection 14.21 is so quantitative and hence we have investigated the difference of the estimated 
value of tnl for the observation of XMMU J2235. 3-2557 between the case with MVJ expression and that with Eq. (jSTJ given in section[3] 
Our naive estimated value of tnl given in subsection I4.2l is tnl = 1-7 x 10®. For the case by making use of MVJ expression, we obtained 
Tnl = 1.1 X 10®. These values seem to be same order and hence the result does not extremely change. 



